Abstract. We have revisited the cahbration of the visual 
binary system CHercuhs with the goal to give the seis- 
mo logical properties of the GO IV sub-giant ^Her A. The 
sum of masses and the mass fraction are derived from the 
most recent astrometric data mostly based on the HIP- 
PARCOS ones. We have derived the effective temperatures, 
the luminosities and the metallicities from available spec- 
troscopic data and tycho photometric data and calibra- 
tions. For the calculations of evolutionary models we have 
used updated physics and the most recent physical data. 
A minimization is performed to approach the most re- 
' liable modeling parameters which reproduce the observa- 
p\| , tions within their error bars. For the age of the ( Her bi- 
^ nary system we have obtained t^Hcr = 3 387Myr, for the 
^ . masses to^^hoiA — 1.45 M0 and m^HcrB = 0.98 A/0, for 
C/^ ■ the initial helium mass fraction Yi = 0.243, for the initial 
Qv^ , mass ratio of heavy elements to hydrogen (f )i = 0.0269 
psj ' and for the mixing-length parameters A^HerA = 0.92 and 



o 



HerB = 0.90 usiug the Canuto & MaziteUi ( |1991| , |1992D 
convection theory. Our results do not exclude that ^ Her A 
^ , is itself a binary sub-system has been suspected many 
' times in the past century; the mass of the hypothetical 
, unseen companion would be m^Hera ^ 0.05 Mq, a value 
' significantly smaller than previous determinations. A cal- 
ibration made with an overshoot of the convective core of 
CHerA leads to similar results but with a slight increase 
of w -f250Myr for the age. The adiabatic oscillation spec- 
trum of C Her A is found to be a complicated superposition 
of acoustic and gravity modes. Some of these waves have a 



Oh 



Q dual character. This greatly complicates the classification 
^ , of the non-radial modes. For £ = 1 the modes all have both 
^ ' energy in the core and in the envelope; they are mixed 
modes. For £ = 2,3 there is a succession of modes with en- 
ergy either in the core or in the envelope with a few mixed 
, modes. The echelle diagram used by the observers to ex- 
tract the frequencies will work for £ — 0, 2, 3. The large 
5^ , difference is found to be of the order of Ai^q « 42 /iHz, in 



agreement with the Martic et al. ( 2001 ) seismic observa- 
tions. 
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1. Introduction 

C Herculis (40 Her; BD +31 2884; STF 2084; 
ADS 10157; IDS 16 375+3147; WDS 16413+3136; 
HR6212; HD 150680; HIP 81693; a = 16^41^17% 
S = +31°36'10" (2000)) is a well known bright visual 
and single lined spectroscopic binary system of naked-eye 
brightness. According to the CDS Simbad data base 
the system is composed of a 2.90 V magnitude GOIV 
sub-giant star and by a 5.53 V magnitude G7V dwarf 
star. The binarity was discovered by Herschell as early 



Fifth Catalog of Orbits of Visual Binary Stars (Hartkopf 
al. 2000) stipulates: "No evidence in the speckle or 



et 



as 1782 (Aitken 1932) and the system has been carefully 
observed by visual binary observers for more than six 
revolutions back to the first reliable measurements by 
W. Struve in 1826. Several orbital solutions have been 
published from the early nineties to the present day. 
The latest are by Heintz ( |1994| ) and Soderhjelm ( |2000D . 
The orbital elements are well determined. They have 
recently allowed improvements of Sproul Observatory and 
HIPPARCOS trigonometrical parallaxes. About a century 
ago, a duplicity of ^HerA was detected from microm- 
eter and meridian observations (Lewis 1906). A period 
Pab. — 12 yr and a semi-major axis of OAa — 0'.'25 was 
obtained for the sub-binary system ^ Her Aa. Ten years 
later, Comstock (1917) noted that the small irregularities 
of the areal velocity in the orbit have the effect of an 
invisible companion having a period of 18 yr and an 
amplitude less than O'.'l. Later, a careful reanalysis of 
the available astrometric and spectroscopic observational 



material led Herman (1941) to conclude that"/- 7 the 
presence of a third body revolving about the brighter 
component of ( Herculis is not definitely indicated". The 
comparison between the observed position angles and 
distances and their values derived from its astrometric 



orbit allowed Baize (1976) to claim that ^HerA is 
split in two stars of respectively 1.051/© and 0.19 Af© 
with an orbital period PAa — 10.5 yr and a semi-major 
axis OAa = O'.'Oe. McCarthy (|l983|) reported that (Her 
contains at least one unseen companion easily detected at 
2.2 fim by infrared speckle interferometry. A note in the 



Hipparcos data for the large- amplitude third-body orbit 
given by Baize". 



Lebreton et al. (1993) tried to determine both age and 
chemical composition of the system, modeling the two 
components simultaneously, but they did not succeeded 
in modeling the secondary consistently. Based on a precise 
spectroscopic analysis Chmielewski et al. ( |1995| , hereafter 
C95) succeeded in modeling both components. They de- 
rived for the age Her = 4.0+0.4 Gyr and for the masses of 
components respectively 1.3 M© and 0.9 Mq. Since 1995 
the HiPPARCOS's parallax of ^Her has been available; it 
was recently improved by Soderhjelm ( 20001 ). The TYCHO 
& HIPPARCOS magnitudes are also available (Fabricius & 
Makarov 2000); they have been connected to the B, V 



and I magnitudes (Bessell 200C). New improved theoreti- 
cal data are also available, viz. opacities, nuclear reaction 
rates and equation of state. The seismic observations of 
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CHer A recently carried out by Martic et al. (2001) indi- 
cate the presence of solar like oscillations. 

The aim of the present paper is firstly to revisit the cal- 
ibration of the C Her binary system using updated physical 
data and theories and improved astrometrical and photo- 
metrical observational material, and secondly to give seis- 
mological properties of C Her A which will be useful to ex- 
ploit future asteroseismological observations. 

Based on the reasonable hypothesis of a common origin 
for both components, i.e. same initial chemical composi- 
tion and age, the calibration of a binary system consists of 
determining a consistent evolutionary history for the dou- 
ble star, given (1) the positions of the two components in 
the HR diagram, (2) the stellar masses and, if possible, (3) 
the present-day surface chemical compositions. The goal 
is to compute evolutionary models that reproduce the ob- 
servations. The calibration yields estimates for the age, 
the initial helium mass fraction and the initial metallicity 
which are fundamental quantities for our understanding 
of the galactic chemical evolution. Within the error bars 
provided by astrometry, one also determines mass values 
consistent with the stellar structure modeling. According 
to the convection theory applied, one derives values for the 
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Table 1. Relevant orbital elements of the ^Her binary 
system. P is the orbital period, a the semi-major axis, i 
the inclination and e the eccentricity. 



References 



34.45 yr l'.'365 131? 6 0.464 
34.45 yr l'.'33 131? 0.46 



Heintz (1994^ 



Soderhjelm (2000) 



"mixing-length parameter" or "convection parameter" A, 
ratio of the mixing-length to the pressure scale height. 
Once the physics is fixed, the modeling of the two compo- 
nents A & B of a binary system requires a set p of seven 
so-called modeling parameters: 

p = <; t*;TOA,mB,Fi, ( -|- ) ,Aa,Ab 



X 

where t^, is the age of the system, ttia and me are re- 
spectively the masses of components A & B, 1^ is the ini- 
tial helium mass fraction, (■y)j is the initial mass ratio of 
heavy elements to hydrogen, Aa and Ab are the mixing- 
length parameters. There are six observables, namely, the 
effective temperatures Toff a & Tcff b, the luminosities La 
& Lb, the sum of masses S and the mass fraction B (cf. 
Eq. H). Once detailed spectroscopic analyses have been 
performed on the system, the present-day surface metallic- 
ities of stars come as additional observational constraints 
- in the case of ( Her the metallicity is available only for 
the brightest component. 



The paper is divided as follows. In Sect. 2.1, 2.2 and 



2.3, respectively, we collect and discuss the astrometric, 
spectroscopic and photometric observations. In Sect. |^ we 
describe the computation of models and the search for 
the modeling parameters. In Sect. ^, we give the results 
with emphasis on the seismological analysis of CHerA. 
Section ^ is devoted to a discussion and finally, we sum- 
marize our results and conclude in Sect, n 



2. Observations of the visual binary C, Her 
2.1. Astrometric data 

The knowledge of individual masses of each companion 
is one cornerstone of any calibration of a binary system. 
As illustrated by calibrations of a Cen, (e.g. Morel et 



al. 200C), the resulting age is very sensitive to the values 
adopted for the masses of the components. We do not re- 
produce here the large bibliography on relevant astromet- 
ric data (e.g. C95); nevertheless we emphasize the discus- 
sion concerning the estimate of masses. For ^ Her we are in 
the fortunate position of having two precise and indepen- 
dent determinations of the trigonometrical parallax and 
also improved orbital elements; from these data, two esti- 
mates of the sum of masses can be derived independently. 



A first determination is provided by a standard photo- 
graphic parallax (Heintz |1994 ) using an improved relative 
orbit and 152 mid-nights of several decades of Sproul Ob- 
servatory long focus photographic observations. The sec- 
ond one is the recently improved adjustment of parallax 



and orbital elements by Soderhjelm (2000) based on HIP- 
PARCOS data collected during 3.25 yr, old ground based 
observations and recent speckle-interferometry measure- 
ments. Table |l| lists the relevant orbital elements of these 
two recent astrometric orbits; they are so close that we 
can safely adopt the means listed in Table ||. The sum of 
masses is provided by Kepler's third law: 



S — niA + rnB 



ip2 ■ 



(1) 



as usual P (yr) is the period, a (arc sec) is the semi-major 
axis of the relative orbit of the companion with respect 
to the primary, toa and ttib are the masses (solar unit) 
of component A and B respectively and m (arc sec) is the 
parallax. 

The outcome of standard long focus photographic par- 
allax measurements is the so-called relative parallax zurci- 
It must be reduced to tUabsj the absolute parallax, by 
adding a correction representing the dependence weighted 
parallax of reference stars. The correction, of order 1 
to 5 mas, is not accurately known (van de Kanip 
For ^Her the correction is 0'.'0043 ( Lippin cott 
With TUici = 0'.'0974 ± 0'.'0039 (Heintz |l994[ ) that leads 
to a photographic absolute parallax of Wahs = 0'.'102 ± 
0'.'0039 (standard deviation). The improved absolute par- 
allax based on hipparcos data and orbital ground- 
based measurements amounts to Wahs = 0'.'0937 ± O'.'OOOG 



1967 ) 



19811) 



(Soderhjelm 200C). We adopt respectively for the parallax 

and the distance modulus of ^Her: 

tu^^Hor = 0'.'094±0'.'001, dcHcr = -0.134 ± 0.023. 

The sum of masses calculated with Eq. (|^) is: 

S^Hcr = 2.50 ± 0.14 Mq. 

Note that the major source of error on the sum of masses 
remains the uncertainty on zu (Couteau 1978| , 40, VI) even 
with modern estimates of the parallax. For the parallax, 
C95 have retained a value close to the former HIPPARCOS 
result W(^iiei = 0'.'0975; with the same orbital elements, 
this larger value for the parallax leads to a smaller sum of 
masses 5^ Hor = 2.4 Mq. 

C Her is one of the rare binary systems for which the 
mass fraction: 

B=^^, (2) 

TOA + rUB 

can be obtained both by astrometric and spectroscopic 
observational data. A first estimate of B is provided by 
photographic data. As the orbits of the primary and sec- 
ondary are similar with respect to the center of masses, 
the semi-major axis @ of the orbit of the primary is related 
to a by @ = Ba. For the so-called resolved binaries the two 
companions are distinguishable by the detector and @ is 
directly derived from the observations. For unresolved bi- 
naries, the distance between the components is below the 
resolving power of the detector, in which case no separa- 
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tion is possible, and a composite image results. A close 
companion affects the center of light by pulling it toward 
the barycenter. With the photographic technique it is as- 
sumed that the measured position of the blended image 
represents the weighted center of light-intensity, or pho- 
tocenter, of the components. In this case, the fractional 
distance /3 of the primary to the photocenter, in terms 
of the distance between the two components, is given by 



(vandeKamp 1967): 



1 



lAAr. 



La + Lb 1 + IQO 
Here La and Lb are the luminosities of the components 
and Am the difference in magnitude. The fractional differ- 
ence of photocenter to barycenter is therefore B — (3 and 
the semi-major axis of the photocenter orbit relative to 
the barycenter is {B~/3)a. This /3 correction also explicitly 
affects the derivation of the parallax. It has long been sus- 
pected that for a magnitude difference Am ^1.5 the posi- 
tion of blended images generally cannot be considered by 
the simple geometric Eq, 



(Morel |197C| , Feierman |l97lD . 
Likely, /3 is also a function of the separation on the photo- 
graphic plate, among other parameters. For the determi- 
nation of the mass fraction of CHer with Am = 2.6 (Hof- 
fleit & Jaschek p82|) Eq. (P) lead to /? = 0.084 while, from 



the work of Feierman ( |1971| ), /? = 0.025 (Lippincott |198lD . 
These differences in (3 lead to the discrepancies between 
the mass fraction derived by Heintz, B — 0.412 ± 0.020, 
and Lippincott, B — 0.352 ± 0.020, despite almost identi- 
cal B — /3 valu es, namely 0.328 & 0.32. We shall refer to 
Heintz's ( |l994| ) value of Bp = 0.41 ± 0.02 as the "photo- 
graphic" mass ratio. 

The second estimate of the mass fraction is provided by 
the spectroscopic orbit. For a single lined spectroscopic bi- 
nary, with known orbital elements and parallax, the mass 
fraction is given by (e.g. Heintz 

B ' 



KawPVT^ 



1971 



Scarfe et al. 1983) 



(4) 



6.283 asini 
As usual i is the inclination, e is the eccentricity and Ka, 

in AUyr^^, is the velocity semi- amplitude of the primary 

star. Combined spectroscopic observations covering about 

four decades and astrometric measurements give Ka — 

4.01±0.04kms^^ (Scarfe et al. |1983D . From spectroscopy 

and orbital elements we obtain the "spectroscopic" mass 

fraction: 

Bs = 0.38 ±0.02. 

We note the good agreement between spectroscopic and 
photographic determinations. 

The mass f ractio n is not derived in the analysis of 
Soderhjelm ( 2000 ) despite an angular distance p ^ 1'.'50 



lar ger th an the grid step (1'.'208) of hipparcos (Martin et 
al. |l997|) , and changes of A9 - 10° and Ap - O'.'IO respec- 
tively in position angle and separation along the flight of 
the satellite. 

For further investigations we shall adopt as the mass frac- 
tion of ( Her the weighted mean: 
B(;Hor = 0.40 ±0.02, 



between astrometric and spectroscopic values. C95 have 
adopted a slightly larger value Bquct = 0.42. 
With the values retained for 5^ Her and B(^ Her the individ- 
ual masses are: 

m^nerA = 1.50 ±0.16, m^HerB = 1.00 ±0.08. 

Table || lists the astrometric data and the astrometric con- 
straints we use for the calibration of the C Her binary sys- 
tem. 



(3) 2.2. Spectroscopic data 



Effective temperatures and luminosities can be deter- 
mined from both photometric and spectroscopic analyses. 
As a general rule, when available, it is the detailed spec- 
troscopic analysis, using hydrogen line profiles, which pro- 
vides the best estimates for effective temperatures, while 
bolometric magnitudes and luminosities are more accu- 
rately determined using photometric data and calibra- 
tions. 

With classical 1.5m to 2.0m telescopes CHer appears 
as a single star under average seeing conditions because 
of the large magnitude difference and of the small an- 
gular distance. Therefore isolated spectra of each com- 
ponent cannot be obtained. The available spectroscopic 
data in C95 allows the derivation of effective tempera- 
ture only for CHer A. We start with an estimate of Ted a 



from Magain's (1987) calibration, assuming a metallic- 
ity of [Fe/H]=0.0 and color index (B-V)=0.65 (see Ta- 
ble 1^). Then we perform a standard lte detailed anal- 
ysis using the curve of growth technique with models 
atmosphere from Gustafsson et al. (1975). Equivalent 



widths used are from C95 and the oscillator strengths 



are from Thevenin (1990). The microturbulence is fixed 
to 1.5kms~^. With the updated solar ir on ab undance 
logepo = 7.46 (Holweger |l979| , Asplund |2000|) we de- 



rived the same metallicity [Fe/H]=:0.04 dex as C95. An 
ionization equilibrium is obtained from curves of growth 
of Fel and Fell which permits us to derive the surface 
gravit y log g = 3.85. We correct the [Fe/H] value in Mag- 
ain's ( 1987| ) formula and derived an improved temperature 
value of Tcff A — 5 820 ± 50K and re- iterate the curve of 
growth analysis. Finally we deduce loggA = 3.75 ± 0.15 
and [Fe/H]A = +0.04 ± 0.03. With these improved in- 
put parameters the scattering of the curve of growth de- 
creases and is satisfactory. C95 have derived a slightly 
smaller value for the gravity log (/a = 3.65 ± 0.20 invoking 
non-LTE effects. Such stars with solar abundances are not 
suspected to suffer from nlte overionization (Thevenin & 
Idiart 199S| ), therefore the surface gravity of CHer A can 
be considered as being well determined. It is remarkable 
that this new effective temperature value we derived, us- 
ing both tycho's photometry (Fabricius & Makarov 2000] ) 
and spectroscopic analysis, is very close to 5 825 ± 40 K, 
the value obtained by C95 using both the photometry 
available in 1995 and the profile of Ha corrected by the 
presence of the light of the secondary. The uncertainty 
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on [Fc/HJa has been estimated by varying Teff and logg 
in their extremes and with the accuracy of the fit re- 
sulting from the quality of equivalent widths. The pre- 
cision obtained for the metallicity is compatible with a 
value of 0.05 dex, which is the more pessimistic estimate 
in C95. The precision obtained for the effective tempera- 



ture results from the partial derivative of Magain's (1987) 
temperature formula combined with the abundance un- 
certainty and errors of the tycho magnitudes listed Ta- 
ble ||. For the effective temperature of CHer A, C95 have 
obtained ±40 K, a similar estimate. 

The projected rotational velocity vsini of CHerA 
amounts to 3.9kms~^ (Fekel 1997). Assuming a parallel 
axis for rotation and orbital motion, sini « 0.655, CHer A 
is therefore a slow rotator. We can safely infer that it is 
the same for the less massive B component. So we can ne- 
glect the small rotational velocity in modeling the internal 
structure of both components. 



Table 2. Photometric data from TYCHO and HIPPAR- 
COS and derived bolometric corrections, magnitudes and 
luminosities. 





CHerA 


CHerB 


Bt 


3.70 ±0.01 


6.39 ±0.01 


Vt 


2.98 ±0.01 


5.47 ±0.01 


Hp 


3.022 ±0.002 


5.705 ± 0.027 


B 


3.539 ± 0.055 


6.369 ± 0.081 


V 


2.890 ± 0.033 


5.556 ±0.059 


(B-V) 


0.649 ± 0.088 


0.81 ±0.140 


BCv 


-0.0566 ± 0.009 


-0.168 ±0.026 


Mboi 


2.699 ± 0.065 


5.254 ±0.108 


L/Lq 


6.55 ±0.39 


0.62 ±0.06 



2.3. Photometric data 

We base the photometric analysis on magnitudes and color 
indexes taken from the TYCHO catalogue (Fabricius & 
Makarov 2000D as it provides consistent data for both 
components. For each component tycho's photometric 
system gives magnitudes Vt and Bt with an accuracy 
of the order of 0.01 mag. hipparcos's photometric sys- 
tem gives Hp magnitudes with an accuracy of the order 
of a few milli-magnitude. We use Bessell's (2000) rela- 
tionships between color index (Bt — Vt) and V — Hp, 
to derive the magnitudes in standard Johnson visual 
V and b lue B filters. In view of figures 3 and 4 of 
Bessell's ( ^OOOD paper, the accuracy of calibrations are es- 
timated to ±0.03 m ag. Th e bolometric corrections, derived 
from Bessell et al. ( 1998 ) for each star are respectively: 
BCv(A) = -0.0566 ± 0.008, BCv(B) = -0.168 ± 0.038. 
To be consistent with the detailed spectroscopic analysis 
of C Her A, the effective temperature of C, Her B is derived 
from Magain 's ( p^87| ) formula using the color index (B— V) 
listed Table ^, the metallicity derived for C, Her A and the 
standard gravity of dwarf stars \ogg ~ 4.44. The accuracy 
on TeffB is derived as for Teff a- Resulting from the high 
accuracy of tycho magnitudes, we derive a better preci- 
sion than C95. We obtain ToffB = 5 300 ± 150 K which is 
close to 5 290 ± 300 K, the value obtained by C95. 

Table || lists the bolometric corrections, the bolomet- 
ric magnitudes and luminosities computed with Mboi© = 
4.74 as the solar bolometric magnitude to be used with 
the Besseh et al. ( |1998| ) calibration. The Bessell ( ^000] ) 
photometric calibration allows us to derive the standard 
Johnson B and V magnitudes either from (Bt&Vt) ty- 
CHO magnitudes, or from (Bt&Vt) and Hp hipparcos 
magnitudes. The values obtained from each differ from 
each other by more than is expected from the accuracy 
of measurements. Moreover, fixing the color index (B-V), 
the effective temperatures calculated with different pho- 



Table 3. Adopted observational data and derived con- 
straints for the calibration of Q Her binary system. 



this paper 



Chmielewski et al. (1995) 



p 


34.45 ±0.005 yr 




a 


1'.'35 ± 0'.'02 




i 


131?3±0?3 




e 


0.462 ± 0.002 




Ka 


4.010 ± 0.04 kms"^ 




VJ 


0'.'094 ± O'.'OOl 


0'.'0975 


S( Her 


2.50 ±0.14 Mo 


2.40 Mq 


13( Her 


0.40 ± 0.02 


0.42 


Teff A 


5820 ± 50 K 


5825 ± 40 K 




6.55 ±0.39 


6.194 ±0.285 




0.04 ±0.03 


0.05 ± 0.05 


Toff B 


5300 ± 150 K 


5290 ± 300 K 




0.62 ±0.06 


0.575 ± 0.069 



tometric calibrations (Magain 1987, Alonso et al. 
Flower 



1996 



1996D differ by more than 150 K. This indicates 



that the errors on data derived with the whole procedure 
may not have a standard Gaussian distribution. Therefore 
the uncertainties of quantities derived from photometric 
data are added in absolute value [C^ norm) instead of 
standard quadrature (£^ norm). Table || allows compar- 
isons between the observational constraints retained in 
C95 and in this paper. The main differences are for the 
luminosity values owing to the smaller bolometric correc- 
tions and parallax we have used. 
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3. Evolutionary models 

The mean angular separation between CHerA & B be- 
ing of the order of p — 1'.'5, the distance between the two 
companions amounts to ~ 150 AU. The tidal interactions 
between the two components are then negligible. We can 
safely assume that stars have evolved as single ones, ig- 
noring each other. 

Mo dels h ave been computed using the CESAM code 
(Morel 1997 ). About 600 mass shells describe each model; 
this number increases up to 2100 for the model used in 
seismological analysis. Around 400 and 30 models are 
needed to describe the evolutions of ^Her A and B re- 
spectively. 

Bas ically the physics employed is the same as in Morel 
et al. ( 2000 ). The ordinary assumptions of stellar mod- 
eling are made, i.e. spherical symmetry, no rotation, no 
magnetic field and no mass loss. The evolutions are ini- 
tialized with homogeneous zero-age main-sequence models 
(zAMS) . In the absence of satisfactory treatment of micro- 
scopic diffusion for stars with mass larger than « 1.3 Mq, 
we do not take into account the diffusion of chemical 
species. This important assumption is discussed in Sect. ||. 

3.1. Physical inputs 

The relevant nuclear reaction rates are taken from the 



NACRE compilation (Angulo et al. 1999). We use the ap- 
proximation: 



Fe 
[— ] 



log(^)0 



log 



log 



where {^^) is the iron mass fraction within Z. We use 



the solar mixture of Grevesse & Noels (1993) i.e. (y)g = 
0.024 5. We employed the OPAL equation of state ( Rogers 
et al. 1996 ) and the opacities of Iglesias & Rogers ( 1992 ) 
complemented at low temperatures by Alexander & Fer- 



guson (1994) opacities. In the convection zones the tem- 
perature gradient is computed according to the Canuto 
& Mazitelli (1991, 1992) convection theory. The mixing- 



length is defined as I = AiJp, where Hp is the pressure 
scale height and A is the mixing-length parameter of or- 
der unity. CHerA, more massive than 1.25 Mq, presents 
a convective core along the main-sequence. Following the 
prescriptions of Schaller et al. (1992) we have also cali- 
brated the CHer binary system with an overshooting of 
the convective core of ^HerA over the distance Ov ~ 
0.2 min(ffp, i?co), where i?co is the core radius. The at- 
mosphere is restored using Hopf's law (Mihalas 1978). We 
use the numerical trick of Henyey et al. ( |1965 ) to connect 
consistently the radiative and the convective parts of the 
atmosphere. Hence a smooth connection of the gradients is 
insured between the uppermost layers of the envelope and 
the optically thick convective bottom of the atmosphere. 



It is an important requirement for the calculation of eigen- 
mode frequencies. The radius i?^ of any model is taken at 
the optical depth where T(r*) — Tag, = 0.645 with 
Hopf's law. The mass M^, of the star is defined as the mass 
enclosed in the sphere of radius i?^. The external bound- 
ary is located at the optical depth Text — 10^*, where the 



density is fixed at a standard value in Kurucz's (1991) 
atmosphere models: p(rcxt) = 3.55 10~^gcm^'^. 



3.2. The search of calibration parameters 

The calibration of a binary system is based on the adjust- 
ment of stellar modeling parameters to observational data 
at the age of the system. Fixing the physics, the effective 
temperature Toff, the luminosity L and the surface metal- 
licity of a stellar model have the formal dependences 
with respect to modeling parameters: 

logTcff(*)mod = logToff ( t^;m^,yu, A^ 



H 



log 



Lr. 



log 



L 



Fe 
^H^ 



t^ ; , 



Fe 



Fe 



]i*,A* 



(5) 



^© / mod 

t H \ — ^ ^ 
where the subscript "mod" refers to model values. With- 
out microscopic diffusion almost no change in chemical 
composition due to nuclear reactions occurs in the enve- 
lope of a star as massive as ^HerA so [^]s « and 
Eq. H becomes trivial. The basic idea of the fitting has 
been developed by Lastennet et al. ( 1999| ). To find a set 
of modeling parameters: 

PCHor = JTlAjTOB,!!, ' Aa, Ab| , 

leading to observables as close as possible to the ob- 
servational constraints logTog (C Her A), log ^ ^^^j^" ^'^ ^ , 

logTeff(CHerB), log (MCgiB) ) and [f]3(CHerA) we 

minimize a (i* , "za , tub ,Yi, [ § ] i , Aa , Ab ) functional de- 
fined as: 

s(AUd-[t]s(A)V , 



x~ 



E 

*=A.B 



l0gTcff(*)mod - logTcff(*) 



mod 



fT(l0gTeff(*)) 



(6) 



where the cr's are the uncertainties associated with the 
constraints. For sets of modeling parameters within the 



ranges detailed in Sect. 4.1, we have computed the evo- 
lution of models from homogeneous ZAMS {t = 0) to 
t ^ 4Gyr. Then the 'was computed using Eq. (^) in 
a refined grid obtained by interpolations. We kept for the 
solution the "best" p = p^Her which corresponds to the 
Xmin- After some hand adjustments these best modeling 
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parameters served to compute the models of ^Her A & B. 
Table ^ lists the confidence limits of the modeling param- 
eters of models computed in this paper. The confidence 
limits of each modeling parameter, the other being fixed, 
correspond to the maximum/minimum values it can reach, 
in order that the generated models fit the observable tar- 
gets within their error bars. Owing to the uncertainties 
of the observation material, mainly on the effective tem- 
perature of ^HerB, we do not attempt to improve the 
solution, or to compute more precisely the error bars af- 
fecting the modeling parameters using the singular value 
decom positi on of the design matrix as developed by Brown 
et al. ( |1994[ ). 



4. Results 

4-.1. Calibration of C, Her 

We have computed the evolution of models with initial 
helium mass fraction, initial ratio of heavy elements to hy- 
drogen, masses and mixing length parameters respectively 
in the ranges 0.24 < Y, < 0.27, 0.026 < [f ]i < 0.030, 
1.35 Mq < mcHcrA < 1.65 M©, 0.85 Mq < m^^HcrB < 
1.15M0, 0.9 < A^HerA < 1-1 and 0.9 < A^hcB < LI- 
For the models of ('Her A, as soon as log Toff a < 3.765 the 
evolution is halted. Table ^ lists the set of "best" model- 
ing parameters we adopt for the C, Her A & B models with 
and without overshooting of the convective core of the 
primary. Figure shows the corresponding evolutionary 
tracks in the hr diagram. For the calibration with core 
overshooting of C Her A, we have not recomputed specific 
sets of models to undertake another minimization. We 
fit the observable constraints within error boxes by adjust- 
ments of and (y )i, the other modeling parameters hav- 
ing the values obtained without overshooting. The larger 
value obtained for the age (~ -1-250 Myr) results from the 
larger amount of nuclear fuel available in the overshoot 
convection core. 



Table 4. Calibration parameters of C, Her binary system 
models lying within the uncertainty boxes. The observa- 
tional constraints listed in Table ^ are recalled in parenthe- 
sis. The subscripts 'V and "s" respectively refer to initial 
and surface quantities. 



S(Mq) 
B 
Yi 

A 



age (Myr) 
Teff (K) 
L/Lq 

(#). 
Xi 

Zi 
R/ Rq 
Avq (/iHz) 



CHerA 



CHerB 



1.45 ±0.01 0.98 ±0.02 

2.43 (2.50 ±0.14) 
0.403 (0.40 ± 0.02) 
0.243 ± 0.002 
0.042 (0.04 ± 0.003) 

0.92 ±0.05 0.90 ±0.10 



Calibration without overshoot for ("Her A 

age (Myr) 3 387 ±15 

Toft (K) 5 818 (5 820 ± 50) 5 413 (5 300 ± 150) 

L/Lq 6.72 (6.55 ± 0.39) 0.65 (0.62 ± 0.06) 
(f). 0.0269 ± 0.0005 

Xi ' 0.737 
Zi 0.0198 
R/Rq 2.56 0.915 

Al^(^iHz) 42.2 157. 



Calibration with overshoot for (Her A 
3 625 ± 22 

5 798 (5 820 ± 50) 5 418 (5 300 ± 150) 
6.82 (6.55 ± 0.39) 0.65 (0.62 ± 0.06) 
0.0272 ± 0.005 
0.737 
0.0200 

2.61 0.920 
40.9 156. 



4-.2. Seismological analysis o/CHerA 

The sub-giant star (Her A has a convective envelope 
which may stochastically excite oscillations as is the case 
in the Sun. Some seismic observations of ( Her A have been 
done by Martic et al. ( |200l| ) and seem to indicate a nar- 
row excess of power around a maximum at 675 /xHz in the 
power spectrum. It is thus interesting to consider the seis- 
mic properties of the model selected by the calibration. 

For our (Her A model without overshooting, we have 
computed a set of adiabatic frequencies of the normal 
modes for degrees (. — 0,1,2,3 in the frequency range 
300 to 1 000 /xHz. The model corresponds to an evolved 
star which has burnt all the hydrogen core and presents a 
radiative core and a convective envelope starting at rc ~ 
0.752 i?*, (Mc ~ 0.9929 M*). At age Her = 3 387 Myr the 
convective core which is present during the main sequence 



evolution has disappeared about 440 Myr ago. A zone of 
varying chemical gradient was formed between the outer 
edge of the initial convective core and the center. This /i- 
gradient gives a rapid variation of the sound speed and a 
large value of the maxi mum of Brunt- Vaisala frequency TV 
(e.g. Unno et al. |l989| ): 

g f 1 d\nP d\np^ 



N- _ 

r yFi 91nr d\nr ^ 
of the order of 3 500 /iHz - as usual P is the pressure, 
p the density, Fi the first adiabatic exponent and g the 
local gravity. It acts like a potential well that may trap 
gravity waves. Figure ^ shows the profiles of the Brunt- 
Vaisala frequency N and of Lamb frequencies Li (i = 
1, 2, 3) with respect to the normalized radius r/R^,, viz. 
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Fig. 1. Evolutionary tracks in the H-R diagram for models of ^Her A, & B without overshooting (full), with over- 
shooting of 0.2 Hp for the convective core of C Her A (dot-dash-dot) and with a larger metallicity of 0.05 dex for C Her B 
(dotted) . The open squares correspond to the locus of ( Her B assuming this larger metallicity. Dashed rectangles de- 
limit the uncertainty domains. Top panel: full tracks from ZAMS. The signs denote 1 Gyr time intervals along the 
evolutionary tracks. Bottom left and right panels: enlargements around the observed ^Her A & B loci. 



the corresponding propagation diagram. Recall that the 
Lamb frequency writes: 



eg. Unno et al. ( 1989 ), Vs is the sound velocity. 

Looking at the set of computed frequencies, we see 
that for each degree, the oscillation spectrum is no longer 
composed of two separated sets of modes with acoustic 
(p-modes) and gravity (g-modes) behavior as in solar-like 
stars but it is a complicated superposition of these two 
sets. Some of these waves have a dual character as they 
behave like pressure waves in the envelope of the star and 
gravity waves in the core. This greatly complicates the 
classification of the non-radial modes. 



The distinction between the g- and p-modes can be 
made by considering their normalized integrated kinetic 
energy (or inertia) En/ (e.g. Christensen-Dalsgaard & 
Berthomieu 1991). The energy of the p- modes does not 



depend on the degree. Figure H] plots this quantity as a 
function of the frequency for £ = 0, 1,2,3. The (7-modes 
have a much larger energy than the p-modes. It appears 
that for £ — 1 the modes all have both energy in the core 
and in the envelop, so they are mixed modes. Figure ^ 
upper panel, shows the distribution of kinetic energy den- 
sity of two of these modes along the radius of the star. For 
£ — 2,3 we have a succession of modes with energy either 
in the core or in the envelope with a few mixed modes as 
seen in Fig. ^ lower panel. 
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Fig. 4. Density of kinetic energy pr^Sr.Sr* of two modes with consecutive radial orders for £ — 1 (upper panels) and 
£ — 2 (lower panels) as a function of the normalized radius (p is the density and Sr the displacement). Note that the 
modes £ = 1 are mixed modes while the modes £ — 2 are alternatively g- and p-modes. 



Table 5. Mass (Mq), initial mass fraction of helium, 
initial ratio of heavy elements to hydrogen, age (Myr), ra- 
dius (Rq), effective temperature (K) and large difference 
(/iHz) of CHer A models - without overshooting - com- 
puted with modeling parameters within the error bar of 
Table g 





M 


Yi 




age 


R 


Tcff 




Ag+ 


1.45 


0.243 


0.0269 


3400 


2.58 


5763 


41.6 


Ag- 


1.45 


0.243 


0.0269 


3372 


2.53 


5868 


42.8 


Am+ 


1.46 


0.243 


0.0269 


3315 


2.59 


5811 


41.6 


Am- 


1.44 


0.243 


0.0269 


3415 


2.53 


5808 


42.7 


Ay-f 


1.45 


0.245 


0.0269 


3298 


2.57 


5808 


41.8 


Ay- 


1.45 


0.242 


0.0269 


3332 


2.55 


5811 


42.3 


Az+ 


1.45 


0.243 


0.0274 


3376 


2.55 


5810 


42.4 


Az- 


1.45 


0.243 


0.0264 


3365 


2.57 


5808 


41.8 



The set of frequencies with acoustic behavior can be 
analyzed by classical asymptotic tools. The acoustic char- 
acteristic mean large frequency spacing, also called "large 
difference" , corresponds to a mean of the quasi uniform 
spacing between mode frequencies with the same degree 
and consecutive orders, is given by: 

^=(2/°^) ~42pHz; (7) 
where Rg corresponds to the radius of the outermost shell 



of the model. The scaling of Kjeldsen & Bedding (1995) 
gives a close value: 



39.7 AiHz, 



where m and r are respectively the mass and the radius 
in solar units and Avqq = 134.9 /iHz. The computed fre- 
quencies are fitted to the following asymptotic relation 
(Berthomieu et al. 1993 ): 

Vn,t = vq + Avi{n - no) + at{n - 7io)^. (8) 
With radial order n between 10 and 20 and no = 15 we 
obtain values for the mean large difference at a given de- 
gree listed Table |[ They are close to the large difference 
At'o, cf. Eq. (0), though smaller. These values are to be 
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Table 6. Theoretical global asymptotic quantities (in 
//Hz) describing the p-mode oscillations according to 
Eq. (|) for degrees ^ = 0, 1, 2, 3. 



0.1 



654.0 40.7 0.130 

1 659.5 38.8 0.336 

2 650.4 40.8 0.094 

3 645.0 40.7 0.111 



Fig. 2. Propagation diagram: profiles as functions of 
the nondimensionned stellar radius, of Brunt- Vaisala fre- 
quency N and Lamb frequencies for degrees t — \ 
(dotted), £ = 2 (dot-dash) and £ = 3 (dot-dot-dot- 
dash). The horizontal lines correspond to the two modes 
£ — 2, V — 563.5 /iHz and v — 593.8 /iHz, the first one 
with gravity behavior and the second one with acoustic 
behavior. Figure ^ shows their eigenfunctions. The full 
line indicates the region of the star where the modes are 
propagative. 
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Fig. 5. Frequency differences between acoustic modes of 
same degree and consecutive radial order as function of 
the frequency (in //Hz). 



Fig. 3. Logarithm of the normalized integrated kinetic 
energy E,^j of the modes as a function of frequency (in 
/xHz) for modes of degree £ = (full), £ = 1 (dashed), 
£ = 2 (dot-dash), £ ~ 3 (dotted). Note that all the points 
corresponding to p-modes are on the full line. 



compared to the observations of Martic et al. (2001) who 
derive a value around 42 /iHz from the construction of 
an echelle diagram. Figure || plots the large differences 
Vn,i — Vn-i,e. with respect to frequency. Except for modes 
£ — 1, the points corresponding to high frequency modes 
are close to the same fiat curve around 41 //Hz, with oscil- 
lations due to the rapid variation of the adiabatic expo- 
nent Fi in the helium ionization zone. As in the Sun, the 
(7-modes have small amplitudes at the surface and thus 
they will be hardly observable. Therefore it is probable 



that the echelle diagram used by the observers to extract 
the frequencies, taking account of their asymptotic distri- 
bution, will work for £ = 0, 2, 3 but not for the mixed 
modes £ — \. The fit of the small differences: 

di3{n) = i^n,i=i ~ ^71-11=3. ^ (51^13 + ai(n- uq), 
gives the following results (/tHz): 
Si^02 = 3.963, ao = -0.117, 
5vi3 = 15.820, ai = -1.774. 

Figure ^ shows the asymptotic behavior of 61/02 with 
aligned points and the large dispersion of points for Svis 
due to the mixed character of the modes £ = 1. 

The calibrated model of ( Her A calculated with over- 
shoot of its convective core, has a larger radius inducing 
a smaller value of Ai/q as seen in Table |. 
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Fig. 6. Variation of the small differences do2 (full dots) 
and 1^13 (open dots), as a function of the frequency (in 
fiHz). The full and dotted lines represent their linear fit. 



5. Discussion 

The Achilles' heel of our calibration is the derivation of ef- 
fective temperature and luminosity of ^HerB. It is based 
on the only reliable measurements of B and V magni- 
tudes of each component by tycho and hipparcos. The 
Bessell ( 2000| ) photometric calibration allows us to de- 
rive the standard Johnson B and V magnitudes either 
from (Bt&Vt) tycho magnitudes, or from (Bt&Vt) 
and Hp hipparcos magnitudes. In the case of ^HerB, B 
and V magnitude values obtained in each way differ from 
each other by more than is expected from the accuracy 
of measurements. Moreover, fixing the color index (B-V), 
the effective temperatures calculated with different pho- 
tometri c cal ibrations (Magain 1987 , Alonso et al. 1996| , 
Flower 1996) differ by more than 150 K. Fortunately, the 
scattering is not so large with the bolometric correction 
BCv- As a result, the effective temperature and the lumi- 
nosity of C Her B are hardly known with a relative accuracy 
better than 3% and 10% respectively. As a matter of com- 
parison we have derived the effective temperature and the 
luminosity of ^ Her A with better accuracy, respectively 
0.8% and 6%. 

We are aware that our models suffer from the absence 
of diffusion of chemicals species. We do not introduce it as 
we do not have a satisfactory description of the physical 
process which acts against the too large efficiency of the 
gravitational settling in the envelope of a main sequence 
stellar model with mass larger than M ^ l.SM© i.e. 
without a significant outer convection zone (e.g. Schatz- 
man 1969, Turcotte et al. 1998). The derived metallicity 



of ^ Her A could not be representative of the initial mix- 
ture that formed ^Her A & B. Because of the large differ- 
ence in mass between the components, the diffusion has 
been more efficient in C Her A than in ( Her B and therefore 
the adopted metallicity for ( Her B needs to be increased. 



As the mass of CHcrB is close to the solar one, we can 
safely assume that the change in metallicity has occurred 
at the same rate in the Sun and in CHerB. From ZAMS 
to the age t « 3 400Myr, the metallicity of the Sun is 
increased by about 0.05 dex. Using Magain's ( |1987 ) for- 
mula, an [Fe/Hje increase corresponds to an increase of 
^ +75 K in effective temperature. Figure |^ shows that 
the locus (open square) of C Her B in the hr diagram with 
Tcft B larger by 75 K is closer than before to the evolution- 
ary tracks. To go further one needs a measurement of the 
metallicity of ^HerB. 

With niA = 1.50 M0 and tob = 1.00 Af0 it was not 
possible to obtain simultaneous satisfactory adjustments, 
within the error boxes, for both components. Realistic so- 
lutions are found with tua — 1.45 Mq. That may indi- 
cates that the suspected duphcity of (^HerA is perhaps 
real. In such a case the hypothetical unseen component 
C, Her a will be less massive than previously announced, 
"^CHera ~ 0.05 Mq and will be a brown dwarf or a giant 
planet. Though Table ^ exhibits insignificant small dif- 
ferences between mixing length parameters of ^ Her A & 
B, we obtained values of the order of unity, as expected 



with the Canuto & MaziteUi (1991 1992) convection the 



ory. The differences between the C95 and our calibrations 
mainly result from the difference in distance. In the HR 
diagram the locus of our C Her A model is found in the 
Hertzsprung gap soon after the main sequence, as ex- 
pected for a sub-giant star. In this part of the hr diagram, 
the effective temperature of a star model varies rapidly 
with respect to time. So the age is very sensitive to small 
changes of physics or of modeling parameter. The small 
uncertainty we give for the age corresponds to the cross- 
ing time of the error box in effective temperature. It is 
valid only with the physics we used and with the central 
values of other modeling parameters. In their study C95 
obtained a value for the age: t^Hor = 4.0 ± 0.4 Gyr close 
to ours. We emphasized the fact that the crossing time 
of the C95 error box in effective temperature (±40 K) is 
about twenty times smaller than their estimated accuracy 
in age. 

Table | lists the ages, the radii, the effective tempera- 
tures and the large differences of C Her A models computed 
with modeling parameters within their error bars, but 
for mixing-length parameters. The age of the model Ag— 
{respt. Ag-|-) corresponds to the time elapsed from ZAMS 
to the instant where the effective temperature crosses the 
left {respt. right) limit of the error box in the HR dia- 
gram, namely log Toff = 3.768 {respt. log Tog — 3.761). 
The ages of other models listed in Table || corresponds to 
the time elapsed from ZAMS to the instant where the effec- 
tive temperature crosses the central value log Toff = 3.765. 
The variations of Avq reflect the differences in radius and 
mass of star models. Table ^ shows that one derives un- 
realistic precision for the modeling parameters by defin- 
ing their accuracy in such a way that any combination 
of them, within their error bars, will provide models of 
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CHer A & B within the observational constraints. Owing 
to the weakness of the observing material we do not at- 
tempt to derive the accuracy of the modelling parameters 
using the method developed by Brown et al. ( |1994D . A 
more sensitive modeling parameter is the age due to the 
fast post main-sequence evolution of ^HerA. The accu- 
racies of modeling parameters listed in Table ^ may be 
optimistic. They mean that, for a value of any modeling 
parameter within its interval of accuracy, one can find a 
set of other modeling parameters, each of them within 
its own accuracy limit, in such a way that they will pro- 
vide models of C, Her A & B within the observational con- 
straints. 



6. Conclusion 

Detailed evolutionary calculations of the visual binary 
system C, Herculis have been performed using opal opaci- 
ties and equation of state, NACRE thermonuclear reaction 



rates, Canuto & Mazitelli (1991, 1992) convection theory. 
The effective temperature, surface gravity and metallicity 
of C Her A have been reestimated using published spectro- 
scopic data (Chmielewski et al. 1995, Thevenin [1990 ) and 
new photometric data of the tycho catalogue (Fabricius 
& Makarov pOOOD . The luminosity of (Her A and the 
effective temperature and luminosity of ^ Her B are cal- 
culated according to ty cho and HIPPARCOS photometric 
data and Bessell ( 2000 ) calibrations. The sum of masses 
have been estimated using updated astrometric relative 
orbit and parallax determinations. The mass fraction is 
derived from spectroscopic and astrometric orbits. We 
have determined the most reliable solution within the 
confidence domains of the observable constraints via 
a minimization. Each solution is characterized by 
p = {icHer,mA,TOB,li, [§]i,AA, Ab}, where t^Her is the 
age of the system, mx and tob the masses of components 
A & B respectively, the initial hehum content, 
the initial metallicity and Aa and Ab the convection 
parameter of each star model. We obtained p^^Hor ~ 
{3 387 Myr, 1.45 M©, 0.98 Mq, 0.243, 0.0269, 0.92, 0.90}. 
As a sub-giant, in the HR diagram CHerA is located 
beyond the main-sequence in the "Hertzsprung gap", 
where the star loci move rapidly. So, all other modeling 
parameters being fixed, the locus of C Her A takes only a 
couple of ten Myr to cross over the error box in effective 
temperature. Therefore, fixing the physics^ the modeling 
parameters of the system appear artificially well defined. 
Fixing the masses obtained previously, we have calibrated 
the binary system assuming an overshoot of 0.20 iJp, 
for the convective core of ^HerA. We did not obtain 
significantly different modeling parameters but a larger 
value for the age, namely ~ -f 250Myr. 

The adiabatic oscillation spectrum of ^ Her A is found 
to be a complicated superposition of acoustic and gravity 
modes. Some of these waves have a dual character. This 
greatly complicates the classification of the non-radial 



modes. For 1=1 the modes are mixed modes with both 
energy in the core and in the envelope, they are mixed 
modes. For £ = 2,3 there is a succession of modes with en- 
ergy either in the core or in the envelope with a few mixed 
modes. This will have an implication for the properties of 
the echelle diagram used by observers which will have a 
smooth behavior only for the modes £ = 0,2,3. The large 
difference is found to be of the order of /S.vq « 41 /xHz 
close to the preli mina ry value derived from observations 
by Martic et al. ( |200l| ). 

Six years ago, in the conclusion of their study, 
Chmielewski et al. (1995) requested "[..] to go further, 



(i) a better photometry of component B and (ii) the ex- 
act value of the parallax" . Our work shows that the cal- 
ibration of ^Her remains a difficult task, even with the 
disposal of improved observational data. Among the bi- 
naries to be calibrated with some confidence, ^ Herculis is 
one of the most interesting owing to the difference of evo- 
lutionary state of components resulting from their mass 
difference. Though difficult, CHer is a reachable target for 
modern spectroscopic and photometric apparatus. C Her 
deserves interest to improve the accuracy of the model- 
ing constraints hence, of stellar models and seismological 
analysis. It is all the more desirable because seismic ob- 
servations of ^Her A will hopefully give new constraints. 
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